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We use the cubic complex Ginzburg-Landau equation coupled to a dissipative linear equation 
as a model of lasers with an external frequency-selective feedback. It is known that the feedback 
can stabilize the one-dimensional (ID) self- localized mode. We aim to extend the analysis to 2D 
stripe-shaped and vortex solitons. The radius of the vortices increases linearly with their topological 
charge, m, therefore the flat-stripe soliton may be interpreted as the vortex with m = oo, while 
vortex solitons can be realized as stripes bent into rings. The results for the vortex solitons are 
applicable to a broad class of physical systems. There is a qualitative agreement between our results 
and those recently reported for models with saturable nonlinearity. 

PACS numbers: 42.65.Tg; 42.81. Dp 



I. INTRODUCTION 

The field of spatial pattern formation in nonlinear sys- 
tems has grown significantly in the last decades (see 
reviews [iHSl). In particular, that growth was signif- 
icantly contributed to by the interest in self-localized 
states ("solitons") and their stability in pattern- forming 
systems, both conservative and dissipative ones. 

The motivation of the present work is to achieve a 
quasi-analytical description of the formation of stable 
self-localized structures in spatially-extended lasers. To 
this end, we consider a complex Ginzburg-Landau model 
with the cubic nonlinearity (CGL3), for which an analyt- 
ical chirped-sech localized solution is well known in the 
one-dimensional (ID) setting [13, Ell- While this solution 
is always unstable, it has been shown that an additional, 
linearly coupled, dissipative linear equation can lead to 
its stabilization in coupled-waveguide models, keeping 
the solution in the exact analytical form 0, [l^, . Self- 
localized states in a wide variety of systems described by 
such coupled linear and nonlinear equations, in both ID 
and 2D, was recently discussed in Ref. 
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The physical system which offers a natural real- 
ization of such models is a broad-area vertical-cavity 
surface-emitting laser (VCSEL), coupled to an external 
frequency-selective feedback (FSF). This system has been 
a topic of interest during the last years since it can display 
a variety of localized structures on top of a non-lasing 
background [T5l - [l9j . In this system, the (complex) intra- 
VCSEL field features nonlinear spatiotemporal dynamics 
due to two-way coupling between the optical field and 
the inversion of the electronic population (driven by the 
injection current), while the feedback field obeys a lin- 
ear equation which is linearly coupled to the main equa- 
tion for the intra- VCSEL field. Previous studies have 
modeled the VCSEL-FSF using various approximations 



and producing a rich variety of stable and unstable local- 
ized modes, including the fundamental soliton po| and its 
complex dynamics [2l|, , as well as side- mode solitons 
supported by an external cavity [22j , and vortex solitons 

All the above-mentioned VCSEL-FSF models have 
been numerically investigated including, at various levels 
of the approximation, physically relevant features, such 
as the feedback delay and electron-hole dynamics and dif- 
fusion. On the other hand, the observed phenomena fea- 
ture strong similarities persisting under progressive sim- 
plifications of the model, such as replacing the feedback 
grating with a Lorentzian filter, the adiabatic elimination 
the electron-hole dynamics, and adopting instantaneous, 
rather than delayed, feedback [2^. This observations 
suggest that a simpler underlying model ma y b e intro- 
duced. In this vein, it was implied in Ref. |14| that a 
simplification towards a simple cubic approximation for 
the nonlinearity, could provide such a model of the CGL3 
type. It was also noted that such a cubic approximation 
to the VCSEL-FSF would place it in the same class of 
models as those previously introduced for coupled opti- 
cal waveguides with active (pumped) and passive (lossy) 
cores in Refs. M. JT2I . [isl . [24 . [25| . and for pulsed fiber 
lasers - in Ref. pa ]. 

In this work we demonstrate that such a CGL3 system 
does indeed allow stable and robust fundamental solitons 
in ID, and, which is the basic novel finding, fundamental 
and multi-charge vortex solitons in 2D. We present the 
bifurcation diagram for the fundamental ID solitons in 
our generalized CGL3 model, and establish their stability 
properties. Going over into 2D, we present numerically- 
generated bifurcation diagrams for the fundamental and 
vortex solitons, establish their stability ranges, and ana- 
lyze a relation to the ID solution. To our knowledge, the 
existence of stable 2D solitons and vortices supported by 
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the cubic nonUnearity in the uniform space has not been 
previously reported in any physical context. As concerns 
the stabilization of 2D dissipative solitons by means of 
the feedback, provided by a linearly coupled dissipative 
equation, this approach was first proposed, in terms of 
anisotropic equations of the Kuramoto-Sivashinsky type 
(its 2D modification), in Ref. [23|. The model was sug- 
gested by applications to the flow of viscous fiuid films, 
rather than optics. 



While our primary motivation is provided by the 
VCSEL-FSF, similar soliton phenomena have been found 
in systems such as lasers with saturable _gain and absorp- 
tion f28l - [3^ and with a holding beam [sj-ls^, as well as 
in VCSEL experiments employin g co upled cavities 37 1 
and a built-in saturable absorber [33|, see also Ref. 
for a review. 
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The results reported in this work may have implica- 
tions to both applied and fundamental studies. On the 
one hand, the VCSEL-soliton systems offer a strategy 
for the development of devices for all-optical informa- 
tion processing applications. On the other hand, the re- 
sults constitute a new contribution to the great variety 
of dynamical phenomena described by CGL systems in 
numerous physical contexts. 

The article is organized as follows. In section |TT] we 
present the model previously considered for lasers with 
the FSF and show how it can be reduced to a CGL3 equa- 
tion coupled to a linear one. In section IIIII we consider 
the stability of the zero solution, which serves both as 
the background for self-localized modes and the source 
of pattern-forming instabilities. We then review the dy- 
namics obtained from direct simulations of the CGL3 
reduced model, which reproduces the spontaneous 2D- 
soliton formation, which was reported, in terms of the 
full model, in Ref. [1^, thus justifying the use of the 
reduction to the cubic nonlinearity. In section HVl we dis- 
cuss the analytical ID solution of our CGL3 model and 
produce a bifurcation dia gram similar to those found in 
more complex models 2^, 13, . 



In section |V] we report the most essential new findings 
for 2D solitons. First, we extend the ID analytic solution 
into that for the 2D stripe-soliton family and describe this 
family, with an intention to identify the stripe soliton as a 
limiting case of vortices. Then, using polar coordinates in 
the 2D plane, we find and characterize a family of vortex 
solitons, whose radius increases linearly with the topo- 
logical charge. The stability of the vortex solitons is also 
analyzed in this section, both for the restricted class of 
cylindrically-symmetric perturbations and full azimuthal 
perturbations. While vortex solitons with high values of 
topological charges m are always subject to the azimuthal 
instability (in particular, for m ^ cx) it goes over into 
the longitudinal instability of the stripe) , stable vortices 
with |m| < 3 are found in our CGL3 model. The paper 
is concluded by section IVTl 



II. THE SYSTEM AND MODEL 

Following Ref. [l^l , we start from the model for the de- 
scription of VCSELs coupled to frequency-selective feed- 
back without delay: 



dE 
dt 



dF 
dt 



-k{1 + ia)E + k(1 + ia)^ i+ \ E\'^ 

-iAj_E + F - iixJraE, 
-XF + aXE, 



(1) 



where k is the cavity decay rate, a is the phase-amplitude 
coupling factor, jj is the pump current, normalized to be 
1 at the threshold in the absence of the external feedback, 
is the transverse Laplacian accounting for diffraction 
in the paraxial approximation, ujm is the detuning of the 
maximum of the frequency-selective feedback profile from 
the laser's frequency at the threshold without the feed- 
back, A stands for the width of the frequency filter, and a 
is the feedback strength in imits of k, i.e., the threshold 
is reduced from ^ = 1 at (t = to /i = 1 — (t/k. 

Truncating the Taylor expansion of the saturable non- 
linearity at the third order, Eq. ([T]) is approximated by a 
specific form of the following CGL3 system 



If = goE + g2\E\^E + (d + iD)A^E + F, 



<§ = -XF + aE, 



(2) 



where 



go = k(1 -I- ia){fi - 1) - ioJm, 
g2 = -k(1 + «a)/U, 
a = fiA, 
D = -l, 



(3) 



Note that Kc{go) is the total linear loss (if negative) or 
gain (if positive), and lm{go) plays the role of the effec- 
tive frequency detuning between the laser and the filter 
maximum. Further, Re((72) is the nonlinear loss (if nega- 
tive) or gain (if positive), while Im(g2) represents the self- 
focusing or defocusing nonlinearity. In the general case, 
real parameters D and d account for transverse diffrac- 
tion and diffusion of the field. In the present work we 
mainly consider d = 0, which is relevant to optics mod- 
els in the spatial domain j40l - |43j . but it may be different 
from zero in other physical situations - in particular, in 
the temporal domain [1, [H, [13, [H. 

This approximation of the saturable nonlinearity by 
the cubic expansion is justified by the fact that the 
higher-order nonlinearity, which usually saturates the 
growth of the intensity in lasers without the feedback, 
is no longer the main saturation mechanism in the pres- 
ence of external frequency-selective feedback. Above the 
threshold, the nonlinear term ilm{g2)\E\^ E induces a fre- 
quency shift that, together with the frequency-dependent 
feedback, introduces an effective saturation capable of 
limiting the field amplitude even without nonlinear losses 
[although Re (52) is typically negative for lasers]. 
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Model ([2]) is precisely the CGL3 equation coupled to 
a linear equation. For laser models, there are usually 
specific relations between go and 52; however, in other 
physical situations all parameters of model ^ may be 
independent, providing for the opportunity to study dif- 
ferent behaviors of the solutions. 



III. OVERVIEW OF THE BEHAVIOR OF THE 
SYSTEM 

Linearizing around the zero (non-lasing) solution, the 
evolution of perturbations Se = {SE, 5F) is governed in 
the Fourier space by equation 



-Se{k^) = A/(fci)5e(kx), 
at 



where 



go-(d + iD)kl 1 
a -A 



(4) 



(5) 



The stability of the zero solution against perturbations 
with wavenumber k± is determined by the eigenvalues of 
matrix M. Note that, since the stability of the zero so- 
lution is independent of nonlinearities, the analysis con- 
sidered here is also valid for Eqs. (P) for corresponding 
parameters. 

In Fig. [1] wc show how the marginal stability curve 
changes with detuning a;,„. Taking into account the def- 
inition of go, for d — matrix M depends on uj„i and 
k± only through the combination w,„ -I- Dk^. Since 
D = —1, increasing Wm the marginal stability curve 
translates rigidly to larger values of /c^ as displayed in 
the figure. The mean slope of the instability balloon de- 
pends on the value of a as shown in Fig. [5] An interesting 
property, exploited in previous works, is the existence for 
positive a and for a range of negative values of ujm of a 
region of stability for the zero solution above the off-axis 
emission threshold [See Fig. [Ha) ]. I n this region, stable 
self-localized modes can be found |20l [22j | . In contrast, for 
a < [see Fig. HJa) corresponding to the self-dcfocusing 
case] there is no such region for any value of Wm . The zero 
solution can be stabilized, though, by a nonzero value of 
diffusion d as illustrated by Fig. [31 Note also that both 
for the self-focusing and defocusing nonlinearity the zero 
solution can be stabilized by filtering of spatial Fourier 
modes in the feedback loop, which can be modeled by a 
wavenumber-dependent feedback strength a(k 1 )[4^. 

Different spatiotemporal regimes are possible in model 
depending on values of the parameters. For a set 
of parameters close to those corresponding to Fig. [Ija), 
which are relevant to the dynamics of lasers, the following 
scenario is observed. For the pump currents in the unsta- 
ble region (ju « 0.45), small random perturbations of the 
zero solution grow exponentially, see Fig. |4l[a), leading 
to a complex 2D spatiotemporal pattern, which sets in 
at t « 80 in Fig. [H^a). The instantaneous spatial profile 
of this chaotic pattern is shown in Fig. [SJa) . If, starting 




FIG. 1: (Color online). Shaded are regions of the instability of 
the zero solution against perturbations with transverse wave 
number k± for different values of detuning Um. The black 
line is the marginal-stability boundary, (a) Um ~ —250, (b) 
i^m = 0, (c) iJrn = 250, (d) ojm = 500. Other parameters are 
cr = 60, K = 100, A = 2.71, a = 5, d = 0. 
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FIG. 2: (Color online). The same as in Fig. [T] for different 
values of a: (a) a = —5, (b) a = 0, (c) a = 5. Other 
parameters are a = 60, k = 100, A = 2.71, uJm ~ 500, d = 0. 



from this regime, pump current fi is further increased 
up to a value at which the zero background is stable 
(/i = 0.52) [see Fig. [4ljb)], a transition is observed from 
densely packed filaments to a 2D set of isolated quiescent 
solitons, sec Fig.[5l[b). Initially, the distance between the 
solitons is small, and the interaction among them leads to 
a partial annihilation, see the abrupt fall of the integral 
intensity at i w 350 in Fig.[4l[a). However, when the den- 
sity of solitons becomes small enough, the resulting set 
of quiescent cavity solitons is quasi-stationary, featuring 
very weak interactions. 

Here we have presented the results for uJm = —270, 
for which the marginal stability curve is the same as in 
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FIG. 3: (Color online). The same as in Fig. [T] for different 
values of d. Even in the self-defocusing case (a < 0) , the zero 
solution is stabilized by nonzero diffusion. Here d — 0.0 (a) 
and d — 0.3 (b). Other parameters are a = 60, k = 100, 
A = 2.71, a = -5, o;^ = 250. 
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FIG. 4: (Color online), (a) The total intensity of the 2D laser 
field, following the switch-on at the pump current jj, = 0.45, 
and then the transition to a set of quiescent localized spots, 
achieved by ramping the pump up to /i = 0.52. (b) The 
variation of the pump current in time which gives rise to the 
dynamical picture displayed in panel (a). The parameters are: 
a = 60, K = 100, A = 2.71, a = 5, aj„ = -270, d = 0. 



Fig. [IJa) but displaced 20 units to the left. For values of 
ujm ~ —250 moving solitons instead of quiescent ones are 
observed in the final state. 

A similar scenario is observed in one dimension, which 
is of special interest because exact ID self-localized so- 
lutions are available in the CGL3 system, as we discuss 
below. 



IV. ONE-DIMENSIONAL SOLITONS 

In the case of one transverse dimension (A_l = d'^ /dx"^) 
an exact analytical solution to Eqs. ^ can be found in 



FIG. 5: The instantaneous 2D amplitude profile correspond- 
ing to the regimes presented in Fig. |4l at: (a) t = 200, cor- 
responding to the pump current at which the zero solution 
is unstable and a complex spatiotemporal regime arises; (b) 
t — 400, after having increased the pump current to a value 
for which the zero solution is stable, and the spatiotemporal 
pattern decays into a set of fundamental solitons. 



the form of d [H [13 

- ma,> 



(6) 



Substituting expressions 

^ max 



cosh(A:a;)]~^"*'^e'"*. 

into Eqs. ([2]), we eliminate 
a 



(7) 



and obtain the following quadratic equation for chirp (3 



^Rc\g2(d - iD)] 



(8) 



Im [g2{d~tD)] 

which yields a single physical root, due to the condition 
that the field intensity 



(9) 



lm[g2{d^iD)] 

must be positive. Note that /3 does not depend on linear 
coefficient 50 1 but only on nonlinear coefficient g2 and on 
the parameters of spatial coupling, d and D. Once /3 is 
known, a complex algebraic equation involving uj and K 
is obtained. By separating the real and imaginary parts 
of this equation we obtain 

Rc(go) ctA 1 



(10) 



(11) 



Re (/?) + Re ( /3 

where 

P={l + iPf{d + iD) 

Next, we obtain a cubic equation for lj: 

030;^ -I- 020;^ -I- aiuj + flp = 0, (12) 

where coefficients ao, oi, 02, and 03 depend on the sys- 
tem's parameters and on /?: 



03 = Re(/3), 

a2 = Re(po)Im(/3) - lm(go)Re(,S), 
ai = (a + A2)Re(/3), 
flo = 02^^ + o'AIm(/3). 



(13) 
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FIG. 6: (Color online). The bifurcation diagram for analytical 
solution (j6]). (a) The amplitude (absolute value of the field 
at the center of the fundamental soliton) as a function of 
the pump current, (b) The growth rate of unstable, neutral 
and least-damped stable perturbation eigenmodes versus the 
pump current /i. Point M designates the drift instability of 
the soliton. The parameters are ujm = —250, a = 5, k = 100, 
o- = 60, A = 2.71. 

Equation ([T^ can be solved analytically, but it is more 
practically relevant to solve it numerically, as in Ref. 
[14| . To summarize, the analytical solutions can be con- 
structed according to the following scheme. 

• (i) Solve Eq. ([5]) for /3, and choose the proper root 
to satisfy the positivity of \Emax\'^ as per Eq. ©. 

• (ii) Solve Eq. ([T^ for u! with the coefficients defined 
by Eqs. pTl [T^ . and (3 produced by the previous 
step. 

• (iii) Calculate K using Eqs. ([TUl [TT|) . and w, with 
(3 produced by two previous steps. 

• (iv) Calculate Emax using Eq. 

• (v) Calculate F„iax using Eq. ([7]). 

Once all parameters of solution ^ are determined, the 
stability of this solution can be analyzed following the 
numerical procedure developed in Ref. p3j . 

Fig. [Hfa) shows the bifurcation diagram for the local- 
ized quiescent solitons using parameter values typical 
for the lasers models. Point A corresponds to the pump 
threshold for on-axis emission. For the pump levels /i 
between points A and B the zero background is unstable 
as shown in Fig. [TJa). These points coincide with the 
origin of two branches of localized structures. The two 
branches collide at C and disappear through a saddle- 
node bifurcation. 

Figure H^b) shows the real part of the most relevant 
eigenvalues resulting of the stability analysis of the up- 
per branch [the one connecting points C, M and A in 
Fig. Elja)]. The soliton solution is stable between points 
C and M, and a drift instability appears at point M p^ . 
The instability spectrum is not shown between points A 
and B because the background zero solution is unstable. 
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FIG. 7: Spatiotemporal dynamics of the ID soliton in the 
course of the development of the drifting instability. The 
parameters are as in Fig. [S] for jj, — 0.7. 

The lower branch of soliton solution connecting points B 
and C in Fig.[6][a) is entirely unstable, as usual [23. [25|. 

Direct simulations starting from the ID analytic soli- 
ton in the unstable region (to the left of point M in Fig.lH]) 
shows the development of the drift instability, sec Fig. [T] 
Notice that once the drift instability sets in the soliton 
moves away from its original location at a constant speed. 

The overall scenario is very similar to that found in 
previous numerical works for the saturable nonlinearity 
|20l [2^ [23j . suggesting that the present CGL3 system 
indeed represents a simple underlying model which cap- 
tures the essential features of more realistic, but also 
more involved, models. 



V. TWO-DIMENSIONAL SELF-LOCALIZED 
SOLUTIONS: STRIPES, FUNDAMENTAL, AND 
VORTEX SOLITONS 

In 2D (A_L = ^ + the ID solution ^ can be 

generalized to a continuous family of the stripe-soliton 
solutions, parameterized by transverse wavenumber kyi 

E = S„,ax [cosh(i^a;)] e''^t^^k,y^ 

F = J;„ax [cOsh(Xx)]"^"''^ e.'^tf.^kyy^ \ > 

The only difference from the above ID solution is a mod- 
ification of linear coefficient 50 j which is replaced by 
go = gQ — {d+iD)ky. This solution is shown in Fig.[8l(a,b). 
We have found the whole family of the stripe solitons as 
a function of ky, see Fig. [5] They exist only for values 
of ky below a certain value beyond which the frequency 
shift introduced by ky pushes the solution outside the 
frequency range of the feedback filter. The solution with 
largest amplitude, marked by a filled circle, corresponds 
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FIG. 8: The spatial profile of the amplitude (a) and real part 
of the field (b) for the unstable 2D stripe soliton (|14p . (c) and 
(d): The same as (a) and (b) for a stable vortex with m = 3. 
Here fi — 0.52 and other parameters are the same as in Fig.[4l 



to ky such that w 



0. Fig. |ni[c), shows the largest real 
parts of the eigenvalues obtained from the linear stability 
analysis in the x-direction. The stripe-soliton undergoes 
a drift instability similar to the one of the ID soliton de- 
scribed in the previous section. Here the stripe as a whole 
would start to move either to the left or to the right of 
its axis. Interestingly enough the drift instability takes 
place at a value of ky = ky which, within the numerical 
accuracy, coincides with the value for which the solution 
has w = 0. The critical value ky will be very relevant 
later when studying the radial dynamics of vortices. In 
any case, 2D stripe-solitons are always unstable to per- 
turbations in the y-direction, breaking up into a number 
of fundamental (spot) solitons. 

We now proceed to fully localized 2D solutions. There 
are two types of stable 2D modes: fundamental solitons 
with the bell-like intensity profile, see Fig. ^h), and 
ring-shaped vortex solitons, see Fig.[5l^c). Vortex solitons 
with integer topological charge m can be looked for as 
E{r,(l)) = Eo{r)e^'^'^, where (r, </>) are polar coordinates 
with the origin at the pivot of the vortex [13, iM, .45. - .48i | . 
The fundamental 2D soliton corresponds to m = 0, with 
the maximum at the origin. Every m vortex has a mirror- 
image — m vortex, therefore for the sake of simplicity in 
what follows we will consider vortex solitons with m > 0. 

Figure [8] shows the similarity between a vortex mode 
and the stripe soliton. Roughly speaking, the vortex may 
be considered as a stripe bent into a closed circle, at least 
for large values of m. Following this similarity, we study 
the radial dynamics of vortices using the radial version 
of Eq. (0), with 
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FIG. 9: (Color online). The continuous family of 2D stripe 
solitons parameterized by ky, for /i — 0.52 and other parame- 
ters taken as in Fig. [l] (a) Frequency uj, (b) the soliton's am- 
plitude |-Emax|, and (c) growth rates of the unstable neutral 
and least damped stable eigenmodes (the stability is consid- 
ered against perturbations depending only on the x coordi- 
nate). 



Using ID analytical solution © to define the initial 
condition as Eini(r) ~ Eo{r = Rq + x), we simulated 
Eq. ([2]) with the Laplacian taken as per Eq. (fTS)) . with 
fixed m. If the initial ring radius Rq is too small, the 
field decays to zero, but for Rq large enough we observed 
the evolution of radius Tmax corresponding to the maxi- 
mum field amplitude towards equilibrium radius R'^^ (m) 
of the vortex soliton of charge m, see Fig. (TU] If Rq 
is much larger than R^^{ni), we observed that the vor- 
tex shrank at a constant speed, which was followed by 
relaxation oscillations as R^^{m) was approached. Ac- 
tually, this is an unusual behavior, very different from 
the typical curvature-driven dynamics [4^, [H^l • The con- 
stant shrinkage speed may be explained by considering 
the quasi-lD dynamics of the stripe-soliton solution. If 
the initial condition has a very large rmax, the ring can be 
considered locally as a stripe-soliton with ky ~ 0, which 
is drift-unstable (see Fig. [9]). The overall curvature of 
the ring breaks the left-right symmetry the stripe had 
in the direction perpendicular to the axis. The symme- 
try breaking is such that the drift takes place towards the 
center and hence the ring as a whole starts contracting at 
a constant speed. The first stage of the dynamics shown 
in Fig. [TO] (at t < 7) is, thus, essentially the same as in 
Fig. [T] As radius rmax shrinks, the effective wavenumber 
ky increases, eventually suppressing the drift instability 
and the ring relaxes to the stable radius R^*^{m). 

Running the simulations for different (large enough) 
values of to, we have produced the dependence of the 
equilibrium radius R^*' on the topological charge to, which 
turns out to be linear, see Fig. [TT] Therefore, vortex 
rings expand as m increases, tending towards the stripe 
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FIG. 10: (a) The spatiotemporal dynamics of the radial profile 
of the field amplitude relaxing towards the equilibrium radius 
of the vortex with m = 2. Here /i = 0.52 and other parameters 
are as in Fig. |4] (b) The dynamics of position rmax of the 
maximum of the field. 

solution in the limit of m — > oo. The inverse of the slope 
of the line in Fig. [TTl i?*^* (to) /m, is the transverse circular 
wavenumber kc, which is, evidently, nearly constant for 
all vortices. The value of kc turns out be very similar 
to the drift-instability critical wavenumber of the soliton 
stripe, ky considered above. 

The mechanism leading to 2D stable vortices discussed 
here has no counterpart in simple curvature driven dy- 
namics of fronts connecting two equivalent states (49l.[5G|. 
In these systems, ID fronts in a 2D system may be sub- 
ject to modulational instabilities but not to drift ones. 
Therefore there is no transient regime in which the ring 
radius changes at a constant rate. The existence of the 
ID soliton drift instability plays a critical role in the dy- 
namics of 2D solitons and determines its stationary size. 

The radial equation allows one to study the radial dy- 
namics independently of the presence of azimuthal in- 
stabilities. In fact, as in the case of the stripe soliton, 
vortices with large to are azimuthally unstable in the 
full 2D problem. The curvature can, however, prevent 
the azimuthal instability for small topological charges, 
and vortices may be stable up to a certain value of to 
[i^li^. Figure El c,d) shows, for instance, a stable vor- 
tex for TO = 3. 

Finally, following the method described in Ref. [23| . 
we have computed the bifurcation diagrams of the soli- 
tons with m = 0, 1 and analyzed their stability, see Fig. 
[T2] The fundamental soliton (vortex) is stable between 
points Mq and Co (Mi and Ci). Point Mo(Mi) again 
corresponds to the onset of the drifting instability of the 
state as a whole. The branches connecting points B and 
Co {B and Ci ) correspond to the solitons which play the 
role of the unstable separatrix. The fundamental 2D soli- 




m 



FIG. 11: (Color online). The equilibrium radius R"*" of the 
vortex solitons versus the topological charge. Parameters are 
as in Fig. 1101 Filled circles mark even integer values of m, 
while the solid line is a linear fitting. Figure [TOl corresponds 
to the leftmost point in this figure. 

ton and the to = 1 vortex have quite a similar bifurcation 
diagram. As compared to ID solitons, see Fig. [51 points 
A and B correspond to the limits of the region where 
the background zero solution is unstable to homogeneous 
perturbations and therefore are the same, however here 
points C and M are located at a lower pump value. 

In addition, we have observed stable vortices with 
TO = 2 and 3. The region of their existence is almost 
identical to the existence region of to = 1 vortex, while 
the stability region is narrower and lies inside MiCi in- 
terval of Fig. [H 

In the case of saturable nonlinearity, system ([1]) gives 
also rise to fundamental 2D solitons and vortices as en- 
countered here [l^l ■ The region of existence and the bi- 
furcation diagrams are quite similar to the ones shown 
here. This is a clear indication that the present CGL3 
model is indeed relevant for other systems, for which the 
analysis in terms of exact stripe solutions is not possible. 

VI. SUMMARY 

In this work, we have introduced the system of cou- 
pled cubic complex Ginzburg-Landau equation and ad- 
ditional dissipative linear equation as the model of laser 
cavities with the external frequency-selective feedback. 
We have observed a qualitative agreement with the re- 
sults recently obtained in models with the saturable non- 
linearity [20., 22,, 23^. In particular, the stability of the 
fundamental 2D solitons is obvious in the case of the sat- 
urable nonlinearity, while it is a nontrivial finding in the 
cubic model. Using analytical considerations and numer- 
ical analysis, we have shown that 2D vortex solitons can 
be interpreted as stripe solitons bent into rings (as il- 
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FIG. 12: (Color online), (a) The amplitude of the soliton with 
m = (dash-dotted line) and m = 1 (solid line) versus the 
pump current, (b) The growth rate of the 6 localized modes 
with highest A, including unstable, neutral and least damped 
stable eigenmodes for the solution branch CiMiA in panel 
(a), (c) The same as (b) but for the solution branch CqMoA 
in panel (a). Parameters are the same as used above in this 
section. 



lustrated by Fig. [8]). This correspondence is clear for 
m oo, but it actually holds too for rather small to, 
since the circular transverse wavenumber appears to be 
the same for all m (see Fig. [TT|) . In such a way, we have 
established the connection between ID and 2D solitons. 



In our system of the coupled cubic and linear equa- 
tions, we have found stable ID and 2D solitons, including 
vortices. These results are important, as they show that 
this system may be considered as the fundamental model 
underlying a wide class of stable soliton lasers. The sim- 
plicity and flexibility of the linear coupling has previously 
been shown to provide the existence of stable ID solitons 
in the models of dual-core waveguides d, [H, [H, [H, HI] . 
Our extension of this approach into the spatial domain, 
and into 2D, means that models of this class may be use- 
ful to describe and analyze stable self-localization for a 
wide variety of physical systems. 
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